In [4]:
#tutorial lists
#https://techoverflow.net/2021/04/25/how-to-make-cartopy-coastline-or-border-lines-thicker-line-width/
#https://matplotlib.org/stable/api/_as_gen/matplotlib.axes.Axes.streamplot.html#matplotlib.axes.Axes.streamplot
rect_white = mpatches.Rectangle((X_coord.min(), Y_coord.min()), np.round(np.ptp(X_coord[0]),2), np.round(np.ptp(Y_coord[:,0]),2),
ec="none", fc='xkcd:light grey', alpha=0.8)
bulan = ['Jan','Feb','Mar','Apr','Mei','Jun','Juli','Agu','Sep','Okt','Nov','Des']
mountains = mpatches.Polygon(mountains_bnd, fc='xkcd:light brown', zorder=1, alpha=0.5)
bgmap = GoogleTiles(desired_tile_form='RGB', style='terrain', url='http://mt0.google.com/vt/lyrs=p&hl=en&x={x}&y={y}&z={z}')
#bgmap = Stamen('terrain')
#plt.rcParams['axes.linewidth'] =3
extent_customiser = [118.6, 120.6, -1.8, 0.2]
fig = plt.figure(figsize=(10,10), dpi=65)
ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
ax.set_facecolor('xkcd:light blue')
#ax.set_extent([X_coord.min(), X_coord.max(), Y_coord.min(), Y_coord.max()], ccrs.PlateCarree())
ax.set_extent(extent_customiser, ccrs.PlateCarree())
ax.add_image(bgmap, 8)
coastal_line = cf.NaturalEarthFeature('physical', 'coastline', '10m', edgecolor='black', facecolor='xkcd:lime')
coastal_line_edge = cf.NaturalEarthFeature('physical', 'coastline', '10m', edgecolor='black', facecolor='None')
ax.add_feature(coastal_line, lw=2, zorder=0)
ax.add_patch(mountains)
ax.add_feature(coastal_line_edge, lw=2, zorder=2)
ax.plot(119.8635651874369, -0.9183030171555117, 'o', markersize=20, markeredgecolor='xkcd:red',
markerfacecolor='xkcd:pink', markeredgewidth=2)
ax.text(119.86, -0.97, "Palu", fontfamily='Helvetica', fontweight='bold', fontsize=18, ha='center', va='top', color='xkcd:red')
ax.plot([118.70, 118.85],[-1.70, -1.70], '-', color='xkcd:yellow', linewidth=5, zorder=4)
ax.text(118.775, -1.69, '{:.0f} km'.format(np.round(0.15*111,0)), fontfamily='Helvetica',
fontweight='bold', fontsize=14, ha='center', va='bottom', color='xkcd:yellow', zorder=4)
stream = ax.streamplot(X_coord, np.flipud(Y_coord), np.flipud(u_speed_array[0]), np.flipud(v_speed_array[0]),
density=1, color=np.flipud(speed_mag[0]), arrowstyle='->,head_width=0.5,head_length=0.5',
cmap=cmr.get_sub_cmap('cmr.flamingo',0.10,0.95), transform=ccrs.PlateCarree(), zorder=3,
broken_streamlines=True, linewidth=2)
ax.quiver(X_coord, Y_coord, u_speed_array[0], v_speed_array[0], transform=ccrs.PlateCarree())
#cartopy issue https://github.com/SciTools/cartopy/issues/1738
ax.spines['geo'].set(linewidth=3, zorder=5)
#ax.set_xticks(np.linspace(X_coord.min()+0.2, X_coord.max()-0.2,8), crs=ccrs.PlateCarree())
#ax.set_yticks(np.linspace(Y_coord.min()+0.2, Y_coord.max()-0.2,8), crs=ccrs.PlateCarree())
ax.set_xticks(np.linspace(extent_customiser[0]+0.1, extent_customiser[1]-0.1,8), crs=ccrs.PlateCarree())
ax.set_yticks(np.linspace(extent_customiser[2]+0.1, extent_customiser[3]-0.1,8), crs=ccrs.PlateCarree())
lon_formatter = LongitudeFormatter(number_format='.2f', zero_direction_label=True)
lat_formatter = LatitudeFormatter(number_format='.2f')
ax.xaxis.set_major_formatter(lon_formatter)
ax.yaxis.set_major_formatter(lat_formatter)
plt.setp(ax.get_yticklabels(), rotation=90, horizontalalignment='right', verticalalignment='center')
cbar = fig.colorbar(stream.lines, orientation='horizontal', fraction=.03, pad=0.05, label='kecepatan angin dalam m/detik')
ax.set_xlim(118.6, 120.6)
ax.set_ylim(-1.8, 0.2)
ax.set_aspect('equal')
